Optimal dynamical control of many-body entanglement 
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The preparation of highly entangled many-body systems is one of the central challenges of both 
basic and applied science. The complexity of interparticle interaction and environment coupling 
increases rapidly with the number of to-be-entangled subsystems, rendering the requirements on 
the control of many-body quantum systems ever more restrictive. We propose an approach that 
allows to derive optimal control Hamiltonians in a purely algebraic fashion. These drive a composite 
quantum system rapidly into that highly entangled state which can be created most efficiently for 
a given interaction mechanism, and which bears entanglement that is robust against decoherence. 

PACS numbers: 



Quantum entanglement is a genuine trait of quantum 
mechanics without classical analogue. Whereas its cen- 
tral role as a resource for any task of quantum computa- 
tion and quantum communication is undoubted, its rel- 
evance even in seemingly classical bio-molecules is under 
active, controversial debate [HE]. Quantum correlations, 
i.e. entanglement naturally arises between interacting 
particles. But, at the same time, the coupling of the 
particles to a large number of environmental degrees of 
freedom degrades quantum coherences that give rise to 
non-classical correlations. Thus, the interplay of inter- 
particle coupling and decoherence determines the time 
evolution of entanglement, and, generally, limitations to 
the life-time of entanglement due to environment cou- 
pling are ever more restrictive for systems of growing size. 
Therefore, we need tools to actively steer many particle 
quantum systems into highly entangled states, which ide- 
ally should be robust against environmental influences. 
Only with such tools at hand will we be able to pro- 
cess quantum information beyond proof of principle ex- 
periments. Furthermore, this endeavor will improve our 
understanding of entanglement evolution in large, open 
quantum systems, and of its possible role for biological 
functionality [5]. 

In the present paper, we derive coherent control tech- 
niques that allow to exploit inter-particle interactions in 
an optimal fashion and, simultaneously, reduce the in- 
fluence of decoherence, in order to create entanglement 
rapidly, and to slow down its decay. Coherent control, 
i.e. the guidance of a system by an externally applied, 
time dependent control Hamiltonian H c {t), has proven 
very successful for optimal gate-implementation ^ and 
for the suppression of decoherence in noisy environments 
[5J |B], for example by restricting the time evolution to 
decoherence- free subspaces of the Hilbert space [7J, or 
by dynamically decoupling the system from the environ- 
ment |8J. Typically optimal control is based on a target 
functional that is to be maximized by a suitably con- 
structed control Hamiltonian. The generic choice for this 
functional is fidelity, i.e. the success probability of reach- 
ing a predefined target state j9]. However, since entan- 
glement is a non-local property, it is invariant under local 



unitary transformations, so that there is a multitude of 
states with equivalent entanglement properties, that gen- 
erally require different resources for their production and 
have different entanglement decay times [10]. Therefore 
the restriction to specific target states is disadvantageous 
when entanglement is the figure of merit. In this context, 
an entanglement measure as target functional is more fa- 
vorable, and its use very naturally results in the creation 
of the state that is best adapted to the system prop- 
erties - i.e. that is most easily accessible for a given 
interaction, and that has the most robust entanglement 
properties under the specific environment coupling. 

Entanglement measures have already been employed as 
target functionals for two-qubit systems, in the context 
of gate optimization [IT] , or for the optimal choice of an 
initial state to evolve into a strongly entangled state |12j . 
But an extension to higher dimensional, multipartite sys- 
tems is a highly nontrivial step, since the quantification of 
mixed state entanglement is excessively elaborate. Owing 
to the fact that dynamical optimization of entanglement 
implies its repeated evaluation, an absolutely necessary 
prerequisite for the usage of an entanglement measure 
as target functional is the possibility to evaluate it ef- 
ficiently. This requirement is met by an algebraic lower 
bound of the multipartite concurrence c(g) |131I14] , which 
reads 

t(q) = Trg® gA < c(g) 2 , with (1) 
A = 4(P+ - P\ <g) . . . <8> Pf - (1 - 2 1 ~ Ar )P_) , 

where P± are projectors onto the symmetric and anti- 
symmetric subspace of the duplicate Hilbert space H®H, 
and Pi (i — 1, . . . , N) is the projector onto the symmet- 
ric subspace of the duplicate Hilbert space Hi ® Hi of 
subsystem i. Like any entanglement measure, r is invari- 
ant under local unitaries [15] , and it is a valid lower bound 
for the convex roof construction 16J for mixed states — 
i.e. the minimal average entanglement over all pure state 
decompositions of g. This bound is particularly tight for 
weakly mixed states |17j . and exact for pure states, what 
is precisely the class of states that we are aiming at when 
trying to reach high entanglement. 
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The feasibility of any control scheme depends on the 
range of control Hamiltonians that can be experimen- 
tally engineered. Since, in general, one cannot assume 
that interactions can be engineered, we presume that 
control induces only single particle dynamics — i.e. we 
consider local control. Thus, we can phrase our general 
control strategy as follows: which local rotations do we 
need to apply, such that the effect of environment is least 
detrimental — or even beneficial — and such that interac- 
tions are optimal for the creation entanglement? 

A very popular control strategy that underlies algo- 
rithms such as Krotov's algorithm [TH] and GRAPE jl] 
is based on variational methods, in which optimal con- 
trol pulses are reached through many repeated forward 
and backward propagations of the target functional. In 
contrast, our strategy is based on the optimization of 
temporal derivatives of r, and requires only one single 
propagation. The underlying idea is that a maximization 
of f will result in a temporal increase of r. However, f 
is independent of local control fields because of the in- 
variance of r under the induced local unitary dynamics. 
Therefore, we will resort to the maximization of the cur- 
vature t, which in the course of time will result in an 
increase of f , which, in turn, yields the desired increase 
of T. 

Both t and r depend not only on g, but also on its 
time-derivatives 

f = 2Trg><8) qA , t = 2Tr(£<g) g + g) A , (2) 
the first of which is governed by the master equation |19j 
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The second derivative g is obtained through iteration of 
Eq. (|3]). The Hamiltonian in the coherent part of this 
master equation is the sum of the local control Hamil- 
tonian H c , and of the system Hamiltonian H sys , involv- 
ing interactions between subsystems, which are assumed 
to be inaccessible to control. The incoherent part C(g) 
in Eq. (|3| arises from the system's interaction with its 
environment, and typically leads to entanglement decay 

Although the control Hamiltonian H c enters at second 
order via g and via g ® g, the curvature f is actually a 
linear function of H Cl since the quadratic terms in H c 
correspond to local unitary dynamics which leave t un- 
changed. However, linear terms describe the interplay H c 
with the inter-particle interactions H sysi and with the in- 
coherent dynamics C(g). Since this results in global or 
non-unitary dynamics, respectively, f actually does de- 
pend on H c . Physically, this reflects that local control 
induces an evolution towards a state with different dy- 
namical properties. And, even though entanglement can 
not be changed through this local driving alone, the in- 
duced changes in the state's dynamical properties can 




FIG. 1: Four nitrogen vacancy (NV) centers in diamond, cou- 
pled by dipole-dipole interactions of different strengths. Nu- 
merical results obtained in this paper correspond to: Ai^ = 
9.8 MHz, A 3 , 4 = 2.7 MHz, and A 2 , 3 = 1-3 MHz, the largest 
three coupling constants, represented by thick lines. Smaller 
coupling constants are represented by thin lines: Ai,3 = 
0.1 MHz, Ai, 4 = 0.3 MHz, and A 2 , 4 = 0.5 MHz. 

influence the impact of interaction and decoherence on 
its entanglement dynamics. More explicitely, the opti- 
mal control Hamiltonian that maximizes f will steer the 
system towards states for which interactions create en- 
tanglement most rapidly, and for which entanglement is 
most robust against decoherence. 

We can now determine the optimal control Hamilto- 
nian in a straightforward fashion: since f depends on H c 
linearly, this dependence can be phrased in terms of a 
scalar product 

T = J2djT hi + f = X-h + TQ (4) 

i 1 

between the vectors h and X. The former contains the 
tunable parameters of the control Hamiltonian which, in 
the present case of local control, reads H c — J2ij hi+sjC^, 
whith a\ the i-th Pauli matrix acting on the j'-th spin. 
X contains the partial derivatives of r with respect to 
the hi, which can be expressed algebraically in terms of 
the current system state, of the system Hamiltonian, and 
of the Lindbladian C. The additive term tq represents 
the free entanglement dynamics, in absence of control. 
Given a maximal control intensity, the scalar product in 
Eq. Q and therefore the curvature f is maximized by 
choosing h parallel to X. Following this procedure, we 
find the control Hamiltonian that is optimal for the in- 
stantaneous state of a given system in a purely algebraic 
fashion, and its continuous application results in the de- 
sired time-dependent control Hamiltonian that induces 
an optimal entanglement evolution. 

To be specific, we will now demonstrate the efficiency 
of our approach in a system of four nitrogen vacancy 
(NV) centers in diamond [32], as sketched in FIG. [T| 
The ground state electronic spin triplet of an NV center 
is a candidate for a quantum bit, since it can be manipu- 
lated via laser pulses |23j , and read out non-destructively 
through fluorescence measurements |24j . Furthermore, 
its coherence time, reaching up to a few micro seconds 
at room temperature, is only limited by weak coupling 
to the surrounding 13 C nuclear spin bath [25]. All NV- 
centers in a diamond are coupled via the dipole-dipole in- 
teraction H sys — j \,i G \ \- The coupling constants 
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Xi.j decay as 1/r 3 with the distance r between two cen- 
ters. Since NV-centers have fixed positions, the Ajj are 
not tunable, so that we have to resort to the aforemen- 
tioned local control. 

To start with, let us consider purely coherent dynam- 
ics. The time evolution of r is shown in FIG. |2] for a 
separable initial state. If no control is applied (dashed 
line), entanglement oscillates around mediocre values 
on the characteristic time scales of _ff sys (indicated by 
dashed vertical lines), due to dipole-dipole interactions, 
and maximal entanglement is never reached. In contrast, 
after switching on a weak local control, with maximal in- 
tensity \\h\\ = 0.26 Ai, 2, where Ai,2 is the largest coupling- 
constant, the system experiences a considerable increase 
of its average entanglement (dotted line). Yet, oscilla- 
tions still prevail, since the dynamics induced by the con- 
trol Hamiltonian is not fast enough to react to the natu- 
ral system dynamics. Nevertheless, by slightly increasing 
the control Hamiltonian's intensity, we can stabilize the 
amount of entanglement at ever higher average values, 
and an intensity of 1.7 A 1.2 results in the evolution to- 
wards maximal entanglement (solid line) . Oscillations of 
t are flattened out, and once maximal entanglement is 
reached it is essentially maintained. Further augmenting 
\\h\\ yields no substantial acceleration of entanglement 
creation, and saturation sets in. This is due to the fact 
that local control on its own cannot create any entan- 
glement, and the only mechanism that does so is the 
dipole-dipole interaction which determines the minimal 
timcscale on which maximal entanglement can be es- 
tablished. In general, at least n — 1 pairwise couplings 
are needed to create genuine n-partite entanglement, and 
the smallest of these coupling constants determines the 
timescale of £e- In the present case of four spins we have 
= 7r/4A2,3, where A2,3 is the third largest coupling 
constant (FIG. [I). 

Let us now explicitely analyze the dynamics induced by 
dipole-dipole interactions and optimal control according 
to Eq. Q. Since many different states exhibit the same 
amount of entanglement, the entanglement dynamics de- 
picted in Fig. [2] itself does not allow to infer the system's 
time evolution unambiguously. Nevertheless, from the 
observation of the dynamics of many different separable 
initial states we can deduce that, first, spins 1 and 2 are 
entangled by their mutual interaction during an interval 
of t = 7r/4Ai j 2, followed by the formation of a second en- 
tangled pair of spins 3 and 4, after t = tt/AX^^. Genuine 
four-particle entanglement and, hence, maximum entan- 
glement is reached once the spins 2 and 3 — from either 
one of the two already entangled pairs — have interacted 
for a time t — 7r/4A2,3- More explicitly, this stepwise in- 
volvement of an increasing number of spins in the many- 
body entangled state reads 

(|o) + |i)r 4 ^ (|oo) + |ii)) ® (|o) + |i))® 2 ^ 

(|00) + |11))® 2 ^4 |0000) + |llll) + i(|1100) + 10011)), 
where each arrow represents an interaction event, and all 




FIG. 2: Time evolution of the lower bound of concurrence 
r of a separable initial state of 4 spins coupled by dipole- 
dipole coupling as depicted in FIG. [T] The free time evo- 
lution (dashed line), the time evolution under optimal lo- 
cal control with intensity \\h\\ = 2.5 MHz (dotted line), and 
\\h\\ = 17 MHz (solid line) are shown. In addition, the dash- 
dotted line represents the time evolution of r for completely 
unaddressable control pulses, with intensity \\h\\ = 17 MHz. 
Characteristic time-scales of -ff sy s are indicated by dashed ver- 
tical lines. 



states are given up to local unitary transformations. 

Compared to current NMR-techniques, the proposed 
control-scheme is much faster in creating entanglement. 
In |26j . an effective Hamiltonian H c g with uniform cou- 
pling between n spins is implemented in the L = n/2 
collective spin subspace, but in order to avoid leakage to 
other subspaces, the magnitude of the effective Hamilto- 
nian has to be negligible with respect to the energy sep- 
aration from other collective spin states. Therefore, the 
control Hamiltonian has to be chosen such that the effec- 
tive uniform coupling is negligible with respect to the in- 
trinsic dipole-dipole interaction — i.e. \H e s\/\H sys \ = e <C 
1. In the present case of 4 spins, H e g induces an evolu- 
tion into a GHZ-state in icHZ — 327r/eA. In our case, the 
creation of genuine multipartite entanglement is realized 
128/e times faster, since there is no such limitation on H c 
and the full entangling capabilities of H sys are exploited. 
Even the numerically powerful GRAPE-algorithm pro- 
vides, at best, only very modest speed-up with respect to 
the entanglement creation rates achieved by our present 
control strategy: an improvement by a factor of 0.9 is 
achieved for the creation of the symmetrically connected 
4-qubit graph state, while GRAPE yields no speed-up for 
the production of the 4-qubit 2D cluster state [27] . 

The above proves the effectiveness of our control strat- 
egy in the case of coherent dynamics. To test its reliabil- 
ity under more unfavourable circumstances, we want to 
impose additional restrictions. 

First, it may be experimentally challenging to address 
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FIG. 3: Time evolution of the lower bound of concurrence 
r of the system depicted in FIG. [I] for a separable initial 
state under the influence of strong dephasing, at a rate Fi = 
0.02 ns^ 1 . The uncontrolled dynamics is depicted by dashed 
lines, solid lines are obtained by optimal control with intensity 
= 17 MHz. 

all the spins individually, and, instead, a single driving 
field acts on all spins simultaneously. In this case of un- 
addressable control, there are only the three components 
of the control Hamiltonian H c = J^tj hi a l that can be 
adjusted. Nevertheless, control with maximal intensity 
= 1.7 Ai, 2 (same as in the addressable case) results 
in the evolution into a highly entangled state as depicted 
in FIG. [2] by a dash-dotted line. Even though the per- 
formance of these global control pulses is not as good as 
the performance of perfectly addressable pulses, it can 
be improved with higher maximum intensities, and an 
intensity of about \\h\\ — 3Ai2 results in maximum val- 
ues of r (not shown here) that are comparable to those 
obtained when perfect addressability is given. 

Second, unavoidable environmental effects render it 
necessary to consider optimal control of an open sys- 
tem which is subject to decohcrcnce. Here we focus on 
the paradigmatic model of dephasing — i.e. fii — a\ in 
Eq. ([3]). The time evolution of r in a system subject to 
dephasing at rate T = 0.02 /is" 1 is plotted in FIG. [3] 



where the free dynamics in the absence of external con- 
trol (as in Fig. [2| is depicted by a dashed line, and con- 
trolled dynamics is depicted by a solid line. Even in 
presence of this very strong decoherence — attainable co- 
herence times in highly purified diamond are two orders 
of magnitude longer than T _1 |25j — our control scheme 
performs very well, and control pulses with an intensity 
of \\h\\ — 1.7Ai,2 drive the system to almost maximal en- 
tanglement of r — 0.95. Despite the fact that dephasing 
takes its toll on entanglement as time goes by, local con- 
trol stabilizes the amount of entanglement at an average 
value which almost doubles that of the uncontrolled sys- 
tem. This is a consequence of our control strategy to 
maximize f in Eq. ([3| that steers the system in a state 
least sensitive to the Lindblad term in Eq. i.e. most 
robust against environment coupling. 

A further aspect of the performance of control se- 
quences in general is their robustness with respect to ex- 
perimental imperfections such as intensity or phase fluc- 
tuations of the laser fields that are used to implement 
the control Hamiltonian. Whereas a final assessment of 
this feature for our present scheme is beyond the scope of 
the present article, the construction of our control Hamil- 
tonian according to Eq. Q is insensitive to such errors 
at first order. Indeed, tests with, e.g., 10% white noise 
added on top of a control pulse suggest a reduction of 
entanglement by 2% only, in agreement with this expec- 
tation. 

In conclusion, the invariance properties of an en- 
tanglement measure employed as a target functional — 
i.e. its independence of a single constituent's coher- 
ent dynamics — prove to establish significant advantages 
over prior approaches that optimize the target state fi- 
delity Furthermore, by its very construction, our opti- 
mal control of the interplay of the system and environ- 
ment dynamics provides new insight into the static and 
dynamical properties of quantum correlations in large 
and noisy composite quantum systems. Finally, while 
demonstrated here with the help of concurrence, the va- 
lidity of our present approach is by no means limited to 
this entanglement measure: very specific types of entan- 
gled states, as for example required for the execution of a 
certain quantum information tasks |28j can be selected by 
target functionals formulated in terms of suitably chosen 
entanglement measures. 
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